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NEAR-FIELD NOISE COMPUTATION FOR A SUPERSONIC CIRCULAR JET 


Ching Y. Loh and Lennart S. Hultgren 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


Abstract 

A fully expanded, high-Reynolds-number, su- 
personic circular jet of Mach number 1.4 is 
simulated, using a 3-D finite- volume Navier- 
Stokes solver, with emphasis on the near field 
noise. The numerical results are generally 
in good agreement with existing experimental 
findings. 

1 Introduction 

The reduction of aircraft noise is an important techno- 
logical challenge as a result of the environmental limita- 
tions and economical consequences of airport noise pol- 
lution. Jet noise, a major noise component for modern 
commercial aircraft engines, consists of turbulent mix- 
ing noise and, in the case of imperfectly expanded su- 
personic jets, broadband shock-associated noise (as well 
as, discrete screech tones under certain circumstances). 
Jet noise has been a vital part of aeroacoustics since the 
early 1950s and continues to be a topic of many exper- 
imental and theoretical investigations; Refs. [1-3] pro- 
vide a comprehensive discussion and further references. 
A rapid advance in its understanding has occurred during 
the last decade, perhaps in conjunction with the develop- 
ment of computational aeroacoustics (CAA), as well as 
better measurement techniques. The issues, challenges, 
and contributions of CAA are discussed in Refs. [4-6], 
and the Rayleigh-scattering technique [7-9] is an exam- 
ple of advances in experimental methods for high-speed 
jets. 

For supersonic jet mixing noise, it is generally agreed 
that there are two different sources [10-12]. One is as- 
sociated with the convection/evolution of large-scale tur- 
bulence structures in the jet and which produces intense 
sound propagating downstream at angles close to the jet 
axis. Traditionally, this has been referred to as eddy 
Mach wave radiation. The other noise source is asso- 
ciated with the fine scale turbulence in the jet shear flow 
and which produces sound propagating predominantly in 


the lateral (90°) direction [13-16]. The work in this pa- 
per concerns the sound produced by the large-scale struc- 
tures in the jet. 

Substantial progress in numerical computation of jet 
mixing noise associated with the large-scale structures 
has been made during recent years by using large-eddy 
simulation (LES) techniques. Mankbadi et al. [17] per- 
formed an axisymmetric LES computation of a super- 
sonic, high-Reynolds-number jet and determined its far- 
held radiated sound using an acoustic analogy technique. 
Zhao et al. [18], Bogey et al. [19], Uzun et al. [20] com- 
puted jet how and the associated noise by using fully 3- 
D LES at a Mach number of 0.9. The Reynolds num- 
ber in these simulations were low to moderate, however, 
5, 000, 65, 000, and 100, 000, respectively. DeBonis and 
Scott [21] used LES to calculate the hydrodynamics of 
a high-Reynolds-number, supersonic jet of Mach num- 
ber 1.4 [8, 9]. In all of the above work, the numerical 
schemes are of high-order finite-difference type. 

In the present paper, encouraged by the successful 
computation of 2-D axisymmetric and 3-D screech noise 
for round jets [22-25], a second-order (in space and time) 
hnite- volume scheme of low dissipation is adopted for 
LES computation of a supersonic jet of Mach number 
1.4. The scheme is the 3-D CE/SE Navier-Stokes solver 
using a hexahedral grid. 


2 The 3-D Supersonic Jet Noise Problem 

The current computation is for the conditions in the 
Panda and Seasholtz [8, 9] experiments on a circu- 
lar, high-Reynolds-number, supersonic, cold jet at Mach 
number 1.4. In these experiments, a nominally fully ex- 
panded jet is issuing from a convergent-divergent (C-D) 
nozzle. 

In what follows, the exit diameter D of the of the noz- 
zle, the density poo, the speed of sound a^, and the tem- 
perature in the ambient fluid are used as scales to 
make all the dependent variables nondimensional. Us- 
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ing the ideal gas isentropic relations and conservation of 
mass gives that the reservoir and jet exit Mach numbers 
are given by 

M r = (Ae/ArXPe/Pr^+W^Me , (1) 


Mr = 


1 2[(P r /P e )(7-l)/7 - 1] 

( 7 -l)[l-(A e M r )2(P e /P r )2/7] 


( 2 ) 


respectively, where the subscripts ’e’ and ’r’ denote jet- 
exit and reservoir conditions, P is the nondimensional 
static pressure, A r /A e is the area ratio of the nozzle, and 
7 is the constant specific heat ratio. Since the jet is nom- 
inally fully expanded 


Pe = 1 / 7 - ( 3 ) 


Furthermore, it also follows that 


Pe = 1/T e , (4) 

P e = Tr(l + V^)/(l + T^e) , (5) 

«e = M e \fTe , (6) 

Pr = lPr/T r , (7) 

u r = M r yjr~ r , (8) 

where p, T, and u denote the nondimensional density, 
static temperature and axial velocity, respectively. In 
the experiment [8], P r /P e = 3.16, T r = 1.01, and 
A r /A e = 9, which together with (l)-(8) largely define 
the experimental conditions. 


3 Numerical Method 

In this work, the three-dimensional numerical scheme 
is designed to solve an LES version of the compressible 
Navier-Stokes equations. The numerical scheme adopted 
in this work is a variation of the low-dissipation finite- 
volume CE/SE scheme. The CE/SE scheme has formal 
second-order accuracy in both space and time [26, 27] 
and naturally captures shocks and other discontinuities 
in the flow. 


3.1 Computational Domain and Unstructured 
Hexahedral Grid 

The computational domain including the C-D nozzle is 
sketched in Fig. 1. The numerical simulation begins at 
the inlet of the C-D nozzle. The geometrical specifica- 
tion of the C-D nozzle is given in [8, 9]. In nondimen- 
sional form, the C-D nozzle extends from x = —4.9 to 
x = 0 and the external domain ranges from x = —0.2 
to x = 17, where x is the streamwise coordinate. The 
external domain diameter is about 10.8. 

For fully 3-D computations of complicated problems 
such as the current jet noise problem where several dis- 
tinct and disparate length scales need to be resolved [4], 



Figure 1 : Geometry of the computational domain, com- 
prising the interior of the C-D nozzle and an external do- 
main; the nozzle extends into the external domain by two 
internal nozzle exit diameters. 


a large number of computational cells are needed to pro- 
vide sufficient numerical resolution. The choice of an 
unstructured hexahedral grid for the CE/SE method leads 
to a significant reduction of the number of computational 
cells, with an associated reduction in computer memory 
and CPU time requirements, as compared to the com- 
monly used tetrahedral one. The drawback is a modest 
increase in numerical dissipation. 

The unstructured hexahedral grid is generated in the 
following manner [30]: first, a 2-D unstructured quadri- 
lateral grid is generated in a circular region (Fig. 2) and 
then the circular disc is advanced one step in the stream- 
wise direction with step size Ax. The advancing front 
technique is carried until x = 17. The step size Ax 
ranges from 0.02 at the nozzle exit (x = 0) to about 
0.06 at x = 15. Between x = 15 and x = 17 there 
is a buffer zone in which Ax increases exponentially in 
size to provide high numerical damping. Inside the C- 
D nozzle, Ax is much refined. Totally, there are 3.17 
million cells in the computational domain, including the 
the buffer zone. The domain covers the near-field region 
which is the source region for the noise. 

Employing a hexahedral grid instead of a tetrahedral 
one for the current CE/SE scheme helps to reduce the 
number of cells from tens of millions to a few mil- 
lions. Still, the number of computational cells is very 
large and the simulation is implemented using parallel- 
computation techniques. To accomplish this, the un- 
structured hexahedral grid generated for the entire com- 
putational domain is decomposed into subdomains (12 
at present) using the efficient mesh partitioning code 
METIS [28]. Figure 3 illustrates a typical partition of the 
computational domain for the current circular jet noise 
problem. The CE/SE flow solver uses message-passing- 
interface (MPI) library calls to exchanging pertinent in- 
formation between neighboring subdomains during the 
simulation. 
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Figure 2: 2-D quadrilateral mesh in a cross section in a 
3-D hexahedral grid 


3.2 Conservation Form of the 3-D Unsteady 
Compressible Navier-Stokes Equations 

The Navier-Stokes equations for a perfect gas can be 
written in the following nondimensional vector form, 
e.g.[ 29]: 


Ut + F x + G y + H z — 0, (9) 

where x, y, and z are the streamwise, transverse, and 
span wise coordinates, and t is time. The five compo- 
nents of the conservative flow variable vector U are: the 
density, the three momentum components, and the total 
energy per unit volume, i.e., 

Ui = p, U 2 = pu , U 3 = pv , U A = pw, 

U 5 = p/(7 - 1 ) + />(^ 2 + + w 2 )/ 2 , 

where, as before, p, iz, and p denote the density, stream- 
wise velocity, and static pressure, and v and w denote 



Figure 3: Partition of the computational domain for par- 
allel computation 


the transverse and span wise velocities. The flux vectors 
in the x, y, and z directions, F, G, and H , respectively, 
have both inviscid and viscous contributions and the de- 
tails of these flux components may be found in a text- 
book [29] or in a previous paper [30] by the present au- 
thors. Note that in the latter, p is to be interpreted here 
as the nondimensional viscosity. 

The nondimensional equations above could form the 
basis for a direct numerical simulation (DNS). Currently, 
this is not feasible for the high-Reynolds-number flows 
of interest here. The present computations are of the 
LES type, where the effects of the unresolved scales, 
i.e., scales smaller than the grid resolution, are modeled. 
In this case, the equations above are considered the fil- 
tered ones governing the resolved scales, where p and 
p are interpreted as simply averaged and the other de- 
pendent variables are interpreted as Favre, (i.e., density- 
weighted) averaged over the subgrid scale. The simplest 
model for the effects of the unresolved scales, or sub- 
grid scales (SGS), on the motion is obtained by using a 
Boussinesq eddy-diffusivity assumption for the subgrid 
shear stresses and heat flux coupled by Smargorinsky’s 
model for the eddy viscosity (eg. [31]). That is, p occur- 
ing in the viscous contribution to the flux terms in Eq. (9) 
is replaced by 

M = 1 /Re + (CstfiSijSij) 1 ' 2 , (10) 

where Sij = \ is the Favre averaged rate- 

of-strain tensor, 5 = (5x5y5z ) 1 / 3 is a local measure 
of the grid size, C s is a (nondimensional) constant, and 
Re = cloqD / Vqq is the Reynolds number for the compu- 
tation. The inverse Reynolds number 1 / Re = 1.8xl0 -6 
in the present computations, which corresponds to the 
value for the jet Reynolds number, Rej = u e D /v e , of 
1.2 x 10 6 in the experiments [8, 9]. Also, note that the 
simplifying assumption is also made that the subgrid- 
scale Prandtl number equal the laminar one, which is 
taken to have the value 0.72. 
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By considering (x,y,z,t) as coordinates of a four- 
dimensional Euclidean space, £ 4 , and using Gauss’ di- 
vergence theorem, it follows that Eq. (9) is equivalent to 
the following integral conservation law: 

(f I m -dS = 0, m= 1,2,. ..,5, (11) 

Js(v ) 

where S(V) denotes the surface around a volume V in 
E a and I m = (£ m , G m , £ m , C/ m ). These equations are 
discretized and solved using a modified verison of the 
weighted-average CE/SE scheme with a weghting factor 
a — 0 . 8 . 

Note that Fureby and Grinstein [32] point out that 
applying flux limiters to finite-volume methods (as is 
done in the current CE/SE scheme), even in the absence 
of any explicit LES assumptions, effectively leads to 
LES schemes with minimal implicit SGS models. They 
demonstrated through an ‘error’ analysis of a particular 
scheme that the flux-limiters (essentially low-pass fre- 
quency filters) build into the algorithm produces addi- 
tional terms in the equivalent differential forms of the 
momentum and energy equations that can be interpreted 
as the SGS stress tensor and flux, respectively. Hence, 
there is an additional implicit SGS model inherent in the 
particular CE/SE scheme used in addition to the explicit 
assumption in Eq. (10). 

3.3 Initial and Boundary Conditions 

Initially, the flow of the entire computational domain is 
set at the ambient flow conditions, i.e., 

1 

Po = 1 , Po = 

7 

u 0 = 0, v 0 = 0, wo = 0. 

At the exterior upstream inlet boundary, the conser- 
vative flow variables are specified to be the same as the 
ambient flow and their spatial derivatives are held at zero. 
At the nozzle inlet, the reservoir conditions as described 
in §2 are imposed. Of course, v r , w r , and all spatial 
derivatives are held at zero. 

Type I CE/SE non-reflecting boundary conditions 
[33, 34] are imposed at the circumferential and outflow 
boundaries. The no-slip boundary condition is applied 
on all the nozzle walls, external as well as internal. 

3.4 Artificial Forcing 

As is common, e.g. [19], in numerical studies of jet mix- 
ing noise, the simulation allows for a low-level, artificial 
excitation, or turbulence seeding, of the jet shear layer 
near the nozzle exit. This seeding is accomplished by 
adding divergence-free velocity perturbations inside the 
jet shear layer at x = 0.3. The disturbances are tempo- 
rally random and are of the vortex-ring form used in [19], 



where 

^(x,r,0; x p , r p ) = 

16 

exp {—d 2 In 2/ A^) ^ e n cos(nd + <f> n ) (14) 

p n= 1 

and r = (y 2 + z 2 ) 1 / 2 and 6 = sin ~ 1 (z/r) are the ra- 
dial and azimuthal coordiates; x p = 0.3, r p = 0.5 and 
A p = 0.02 define the location and spatial extent of the 
perturbation, and d = \{pc — x p ) 2 + (r — r^) 2 ] 1 / 2 ; e n and 
<fi n are uniformly distributed random numbers between 
— 1 and 1, and 0 and 27 r, respectively. They represent the 
random perturbation amplitude and phase, and are gen- 
erated at every time step. The overall perturbation ampli- 
tude S p is either set to zero (no forcing) or to 0.01 in the 
present computations. The perturbations above are then 
directly added to the velocity field in the latter case. 


4 Numerical Results 

The simulation starts from the uniform initial condi- 
tions described in §3.3, without any artificial forcing. 
The simulation is first run 100,000 time steps, or to 
t — 100 with the time step At = 0.001 used, to en- 
sure that all start-up transients are convected out of the 
computational domain and that the jet flow is fully es- 
tablished. This time, t = 100, corresponds roughly to 
five acoustic times. The acoustic time is a measure of 
how long it takes a sound wave to pass through the ex- 
ternal computational domain. The simulation is then run 
for an additional 340, 000 time steps to t = 440, or ap- 
proximately an additional 17 acoustic times, to achieve 
sufficient accuracy in the spectral analysis of time series 
data that is saved at various locations in the computation 
domain every 20 time steps. The artificial forcing is then 
turned on with an amplitude S p = 0.01 and the simula- 
tion is again run 100, 000 time steps to remove any tran- 
sients. The simulation is then run an additional 350, 000 
time steps, with the random turbulence seeding in the jet 
shear layer in the vicinity of the nozzle lip in effect, and 
time series data is saved at the same locations as in the 
unforced case, again every 20 time steps. The forced and 
unforced results are compared below. 


4.1 Time-Averaged Data 

Figure 4 shows a comparison of the evolution of the 
time-averaged centerline streamwise velocity between 
the simulation (solid and dashed lines) and the experi- 
ments [8, 9]. In general, the agreement is good, particu- 
larly for the forced case, but with the end of the potential 
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Figure 4: Time-averaged centerline stream wise velocity 
evolution; solid line: computation with artificial forcing, 
dashed line: computation with no forcing, symbols: ex- 
periment 

core predicted by the simulation occurring somewhat too 
far downstream. 

Figure 5 shows comparisons of time-averaged radial 
density profiles from the simulations (solid and dashed 
lines) with experimental results [8, 9] for six stream wise 
stations, x = 2,4,. ..,12. The agreement for the first 
two stations (x = 2 and 4, respectively) is very good. 
This indicates that the early jet shear layers are quite well 
resolved in the simulations. The agreement for x = 6 is 
still reasonable, but it is poor for x = 8 and 10, reflect- 
ing the discrepancy also seen in Fig. 4 at these stations. 
However, the agreement is again quite good for x = 12. 
This suggests that the simulation captures the more ener- 
getic/violent mixing process that occurs after the poten- 
tial core is gone. 

4.2 Instantaneous Flow Fields 

Figure 6 shows the computed instantaneous u-velocity 
field in the spanwise plane, y = 0, at t = 440 and 
t = 890 for the no-forcing (top) and random-forcing 
(bottom) cases, respectively. The rather weak shock-cell 
like structures in the upstream part of the potential core 
indicate that the expansion is not 100% perfect as a result 
of either the C-D nozzle design, or the numerical dis- 
cretization/implementation in the nozzle, or both. The 
potential core is observed to end around 10 diameters 
downstream of the nozzle exit and the mixing process 
then rapidly gains strength. Comparing the two plots in 
the figure, the instantaneous flows look qualitatively sim- 
ilar. 

Figure 7 depicts the computed instantaneous spanwise 
vorticity, (j y , in the same spanwise plane, y = 0, and 
times t = 440 and t = 890 as in the previous figure. 





0.90 1.00 1.10 1.20 1.30 1.40 1.50 1.60 1.70 


density 




Figure 5 : Radial variation of time-averaged density at six 
stream wise locations, x = 2,4,. ..,12; solid line: com- 
putation with artificial forcing, dashed line: computation 
with no forcing, symbols: experiment 
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The instantaneous flow in the case with random forc- 
ing (bottom) seems more mixed. Otherwise, the mixing 
processes qualitatively look quite similar in both plots. 

Figure 8 shows the computed instantaneous stream- 
wise vorticity uo x in four streamwise (or cross-sectional) 
planes, x = 3.75, 7.5, 11.25, and 15. The top panel is the 
case without forcing at t = 440 and the bottom panel is 
with random forcing at t = 890. These uo x plots clearly 
show, in both cases, the 3-D nature of the unsteady mo- 
tion at all stations. Again, the unsteady motion is quali- 
tatively similar in both cases. 



Figure 6: Instantaneous streamwise velocity distribution 
in a span wise plane; top, no forcing at t = 440; bottom, 
random forcing at t — 890 


Figure 9 shows the computed instantaneous 3-D pres- 
sure isosurfaces at t = 440 (without forcing) and t = 
890 (random forcing) at the two levels of p = 0.71 and 
p = 0.715. Note that the nondimensional ambient pres- 
sure is 1/1.4. The isosurfaces show downstream radiat- 
ing 3-D sound waves caused by the large-scale coherent 
structures in the jet flow, i.e. eddy-Mach wave radiation. 
In the bottom plot, artifacts of the random forcing can 
be seen in the vicinity of the nozzle exit, but these are 
rapidly washed out in the downstream direction. 

The computations show that the unsteady motion, as 
in the experiments [8, 9], is truly three dimensional. It is 
observed that the flow unsteadiness is sustainable even 
without any explicit random forcing and qualitatively 
similar unsteady structures are still generated. However, 
the discretization and round-off errors that, in fact, sup- 


z 



J — * — ‘ — I — 1 — I — 1 — I — 1 — 1 — I — 1 1 I I I 

5 10 

Figure 7: Instantaneous span wise vorticity distribution 
in a span wise plane; top, no forcing at t = 440; bottom, 


random forcing at t = 890 


ply random forcing in this case are likely of too low an 
intensity to be truly representative of the experimental 
turbulence levels. 

4.3 Spectral Results 

Figure 10 shows the power spectral density (PSD) of 
the density fluctuations at three locations, (x,y,z), in 
the flow field, namely (2,0.48,0), (4,0.41,0), and 
(12,0,0). The first two locations are in the jet shear 
layer, whereas the last one is on the centerline after the 
end of the potential core. The solid and long-dashed lines 
are simulation results with and without artificial forcing, 
respectively. The short-dashed line represents the exper- 
imental results [8, 9]. The Strouhal number, consistent 
with the scaling used here, is defined as S = fD/a oo, 
where f is the dimensional frequency. It is quite com- 
mon to use the jet velocity, rather than the ambient speed 
of sound, in the definition of the Strouhal number — 
dividing our Strouhal number by 1.2 (the present nondi- 
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Figure 8: Instantaneous stream wise vorticity distribution 
in four stream wise planes; top, no forcing at t = 440; 
bottom, random forcing at t — 890 


mensional jet velocity) converts it to the commonly used 
one. 

The top plot in Fig. 10 is for the measuring point in the 
shear-layer closest to the nozzle, (2, 0.48, 0). The broad 
hump in the experimental spectrum represents the hydro- 
dynamic instability waves growing in the jet shear layer. 
It includes large-scale to small-scale disturbances. Due 
to limited grid resolution, and, of course, its general na- 
ture and underlying assumptions, the LES approach can 
not be expected capture the small-scale disturbances, nor 
maybe intermediate-scale ones. This is the reason for 
the sharp drop off in the numerical PSD curves. The 
agreement between the simulation and the experiments 
is, however, quite good for low Strouhal numbers (up to 
about 0.3), i.e., for large scale disturbances. The dis- 
agreement close to S' = 0 is related to the very short 



Figure 9: Instantaneous isobars; top, no forcing at t = 
440; bottom, random forcing at t — 890 


total observation time in the simulation compared to the 
experiment. The simulation results are quite insensitive 
to the artificial forcing. The middle plot in Fig. 10 cor- 
responds to a location in the shear layer further down- 
stream, (4, 0.41, 0). The agreement between simulation 
and experiment is good also here for large-scale distur- 
bances, and again the simulation results do not appear to 
be influenced much by the random forcing. 

The bottom plot in Fig. 10 compares the simulated 
and experimental PSD for the density fluctuations at a 
point on the centerline even further downstream, (12, 
0,0). At this stream wise location, there is no longer a 
potential core and the flow is highly unsteady and three- 
dimensional. While the level of the simulation results 
are too high, the decay rate is comparable to the experi- 
mental data for S less than about 0.6, and then becomes 
steeper. This is an indication that for the current grid and 
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LES approach, the cut-off Strouhal number is about 0.6. 
The much higher level of the numerical results, indicates 
that value of the Smargorinsky’s eddy viscosity constant 
needs to be increased in this flow region, similar to our 
experience in the jet screech computations [30]. Again, 
there is practically no difference between the simulations 
with and without artificial random forcing. 

5 Concluding Remarks 

The present mixing-noise computation for a fully ex- 
panded supersonic circular jet demonstrates that the near- 
field large-scale noise can be determined using a low- 
order (2nd order) finite- volume LES scheme, such as the 
current one. Generally, time-averaged results, such as for 
the streamwise velocity and radial density profiles agree 
well with existing experimental data. Near-field spectral 
results also agree well with the experiments in the ex- 
pected range of lower Strouhal numbers. Instantaneous 
iso- surfaces of pressure show the radiating sound waves 
(eddy Mach wave radiation) originating from the large- 
scale structures in the jet. Instantaneous streamwise- 
velocity and vorticity results demonstrate the highly un- 
steady and three-dimensional nature of the flow. The cur- 
rent results are in good qualitative agreement with what 
is generally understood about the flow physics of high- 
speed jets. 

It is believed that the results can be further improved 
by lowering the damping in the numerical simulation by 
carefully reducing/optimizing the flux limiters and SGS 
parameters in the scheme. Further research in this area is 
currently being carried out. 
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